Spin-plasmons in topological insulator 

D.K. Efimkin\ Yu. E. Lozoviki'^,* and A. A. Sokoliki 

Institute for Spectroscopy, Russian Academy of Sciences, 

Fizicheskaya 5, 142190, Troitsk, Moscow Region, Russia 

^Moscow Institute of Physics and Technology, Institutskii Per. 9, 141700, Dolgoprudny, Moscow Region, Russia 



o 

(N 






CD 



a 



o 



> 

in 

^^ 

\>- 

o 



X 



Collective plasmon excitations in a helical electron liquid on the surface of strong three- 
dimensional topological insulator are considered. The properties and internal structure of these 
excitations are studied. Due to spin-momentum locking in helical liquid on a surface of topological 
insulator, the collective excitations should manifest themselves as coupled charge- and spin-density 
waves. 



I. INTRODUCTION 

In recent years, topological insulators with a non- 
trivial topological order, intrinsic to their band structure, 
were predicted theoretically and observed experimentally 
(see [1] and references therein). Three-dimensional (3D) 
realizations of "strong" topological insulators (such as 
Bi2Se3, Bi2Te3 and Sb2Te3) are insulating in the bulk, 
but have gaplcss topologically protected surface states 
with a number of unusual properties [2]. These states 
obey two-dimensional Dirac equation for massless parti- 
cles, similar to that for electrons in graphene [3], but re- 
lated to real spin of electrons, instead of sublattice pseu- 
dospin in graphene. 

The consequence of that is a spin-momentum locking 
for electrons on the surface of strong 3D topological in- 
sulators, i.e. spin of each electron is always directed in 
the surface plane and perpendicularly to its momentum 
[1, 4]. The surface of topological insulator can be chem- 
ically doped, forming a charged "helical liquid" . 

Collective excitation of electrons in such helical liq- 
uid were considered in [o], where relationships between 
charge and spin responses to electromagnetic field were 
derived. It was shown that charge-density wave in this 
system is accompanied by spin-density wave. Applica- 
tion of spin-plasmons to create "spin accumulator" was 
proposed in [G]. Also the surface plasmon-polaritons un- 
der conditions of magnetoclcctric effect in 3D topological 
insulator were considered [7]. 

In the present article we consider the properties and 
internal structure of spin-plasmons in a helical liquid. 
Within the random-phase approximation, we derive plas- 
mon wave function and calculate amplitudes of charge- 
and spin-density waves in the plasmon state. 



where the Fermi velocity vp w 6.2 x 10^ m/s [-']. Its 
eigenfunctions can be written as e*P'"|/p-y)/-\A5, where S 
is the system area and |/p.^) = (e-''^p/2,J7e''^p/2)T/^ 
is the spinor part of the eigenfunction, corresponding to 
electron with momentum p (its azimuthal angle in the 
xy plane is ipp) from conduction (7 = —1) or valence 
(7 = +1) band. Many-body Hamiltonian of electrons 
populating the surface of topological insulator is -ff = 
Ep7 ?P7ap7«P7 + (1/25') Eq VqP+Pd, whcrc OpT- is the 
destruction operator for electron with momentum p from 
the band 7, ^p^ = 7f f |p| — /i is its energy measured from 
the chemical potential fi, Vq = l-ne^ jeq is the Coulomb 

interaction; p\ = Ep-^7'(./p+q,7'l/p7)«p 
charge density operator for helical liquid. 

The creation operator for spin-plasmon with wave vec 
tor q can be presented in the form: 
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This operator should obey the equation of motion 
[iy, Qq ] = fi^Qq, where $7^ is the plasmon frequency. 
We can get solution of this equation in the random phase 
approximation at T = (similarly to [>^]): 
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where np+ = 6(pF — |p|) and np_ = 1 are occupation 
numbers for electron-doped helical liquid (pp == /i/i'F is 
the Fermi momentum). 

The plasmon frequency is determined in this approach 
from the equation 1 — Vqn(g, fJq) =0, where 
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II. WAVE FUNCTION OF SPIN-PLASMON 

Low-energy effective Hamiltonian of the surface states 
of Bi2Se3 in the representation of spin states {| f), | i)} 
is Hq = VF{px<Jy — PyCr-x) for a surface in the xy plane. 



is the polarization operator of the helical liquid, differ- 
ent from that for graphene [■'>] only by degeneracy factor. 
The factor A^q in (2) can be determined from the nor- 
malization condition 



Qq, Qq' 0) = '^qq' X! ^t'1 "" '^11'' 
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(|0) is the ground state), so that |7Vq|~^ ~ 
— S[dll{q, Lu) / du!)]\uj=nq- The quantities -D-y'-y in (4) can 
be considered as total weights of intraband (£*++) and 

interband {D^ + D ^ = 1 — D++) electron transitions, 

contributing to the plasmon wave function (1). Note 
that all these formulas are also applicable to the case 
of graphenc. 

Spin-plasmon dispersion flq and contribution of intra- 
band transitions into its wave function are plotted in 



Fig. 1 at various 



'/ev-p, where e is the dielectric 



susceptibility of surrounding 3D medium. For Bi2Se3, 
Ts « 0.09 with £ w 40 for dielectric half-space [■'3] (for such 
small Ts, the corresponding dispersion curve approaches 
very closely to the upper bound u = vpq of the intra- 
band continuum). The results for suspended graphenc 
with rather large r^ = 8.8 (for vp « 10^ ni/s, e = 1 and 
with the degeneracy factor 4 incorporated into rg) are 
also presented for comparison. It is seen that the un- 
damped spin-plasmon consists mainly of intraband tran- 
sitions. When the dispersion curve enters the interband 
continuum, the spin plasmon becomes damped and inter- 
and intraband transitions contribute almost equally to its 
wave function. 



III. CHARGE- AND SPIN-DENSITY WAVES 



The helical liquid in the state |lc 



= QJIO) 



with 



one spin-plasmon of wave vector q has a distribution of 
electron-hole excitations (2), shifted towards q. Due to 
the spin-momentum locking, the system acquires a total 
nonzero spin polarization, perpendicular to q. A simi- 
lar situation occurs in the current-carrying state of the 
helical liquid, which turns out to be spin-polarized [4]. 

Introducing one-particle spin operator as s = cr/2, we 
can calculate its average value in the one-plasmon state 

(S) = (lq|s|lq) as 
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FIG. 1: Dispersions of spin-plasmon (a) and contributions 
D++ of intraband transitions into its wave function (b) at 
various rg. Continuums of intraband (a; < vpq) and interband 
{oj + VFQ > 2fj,) single-particle excitations are shaded in (a). 
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If q is parallel to e^;, only the y-component of (s) is 
nonzero. Its dependence on q at various rg is plotted 
in Fig. 2(a). 

Charge- and spin-density waves, accompanying 
spin-plasmon with the wave vector q, can be 
characterized by corresponding spatial harmonics of 
charge- and spin-density operators: pj and s+ = 
Ep77'(/p-i-q,7'|s|/p7)ap+q,y/apy. Using, similarly to [9], 
the unitary transformation, inverse with respect to (1), 
we can write: /9+ = SN^Il*{q,nq)Q'^ + p+ and s+ = 
SN^'n.l{q,nq)Q^ + Sq, where the operators Pq and §+ 
arc the contributions of single-particle excitations and are 
dynamically independent on plasmons. Here the crossed 
spin-density susceptibility of the helical liquid [l] has 



been introduced: 
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The average values of p^ and s+ in the riq-plasmon 

state |nq) = [((5q)"''/(n-q!)~^/^]|0) vanish, therefore we 
consider their mean squares in \nq) after subtracting 
their background values in |0), i.e. 
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= n^S'\N^Uiq,nq)\^, (7) 

(4(4)+) ^ (nq|4(4) + |nq) - (0|4(4)+|0) 

= n^S'\N*nj{q,nq)\^ (8) 

(only the in-plane transverse component s^ of the 
spin s is nonzero in these averages). The normalized 




FIG. 2: Total spin polarization (s) of the helical liquid in 
the one-plasmon state (a) at various rg and normalized am- 
plitudes Ap and As of charge- and spin-density waves respec- 
tively in the many-plasmon state (b). 



amplitudes Ap{q) = [(PqPq)/"^q<5'p]^/^ and As{q) = 
[(*q (Sq)''')/'^q'5'p]"^^^ of cliargB- and spin-density waves 
are plotted in Fig. 2(b) {p = Pp/4:n is the average elec- 
tron density). The "continuity equation" for density and 
transverse spin, following from the spin-momentum lock- 
ing [o], requires that ftqAp^q) ~ 2vFqAa(q), in agreement 
with our results. 



IV. CONCLUSIONS 



We have considered microscopically spin-plasmons in 
helical liquid in the random phase approximation. The 
developed quantum-mechanical formalism can be applied 
for a number of problems in spin-plasmon optics. 

We calculated the average spin polarization, acquired 
by the helical liquid in a spin-plasmon state, as well 
as mean-square amplitudes of charge- and spin-density 
waves, arising in this state. Coupling between these am- 
plitudes, caused by spin-momentum locking, was demon- 
strated. The interconnection between charge- and spin 
density waves can be applied for constructing various 
spin-plasmonic and spintronic devices. 
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